Doping Dependence of Spin Dynamics in Electron-Doped Ba(Fei_2;Coa;)2As2 
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The spin dynamics in single crystal, electron-doped Ba(Fei-a:Coa;)2As2 has been investigated by 
inelastic neutron scattering over the full range from undoped to the overdoped regime. We observe 
damped magnetic fluctuations in the normal state of the optimally doped compound (a; = 0.06) 
that share a remarkable similarity with those in the paramagnetic state of the parent compound 
{x — 0). In the overdoped superconducting compound (a:: = 0.14), magnetic excitations show 
a gap-like behavior, possibly related to a topological change in the hole Fermi surface (Lifshitz 
transition), while the imaginary part of the spin susceptibility x" prominently resembles that of 
the overdoped cuprates. For the heavily overdoped, non-superconducting compound {x = 0.24) the 
magnetic scattering disappears, which could be attributed to the absence of a hole Fermi-surface 
pocket observed by photoemission. 

PACS numbers: 74.70.Xa, 74.72. Ek, 75.40.Gb, 78.70.Nx 



I. INTRODUCTION 



One major difference between conventional and high- 
Tc-cuprate superconductors is the proximity to a com- 
peting magnetically ordered state in the latter, and it 
has long been believed that magnetic fluctuations could 
replace the role of phonons in mediating an electron- 
pairing interaction. This mechanism could give rise to 
more tightly bound Cooper pairs, elevating the transi- 
tion temperature. The recent discovery of iron pnictide 
superconductors^ with Tc exceeding 50 K (Ref. 2) in close 
proximity to antiferromagnetic order reinvigorates this 
idea. 

For the parent compounds of the cuprates, the 
magnetic properties are well described by the two- 
dimensional (2D) quantum non-linear sigma modelf2.ii 
and magnetic order is driven by a large instantaneous 
2D correlation length, weak interlayer coupling, and 
spin anisotropics.^ For the iron pnictides, on the other 
hand, there is still much debate over the nature of 
the magnetism.— For example, it remains controver- 
sial whether the stripe-type antiferromagnetic order in 
the parent compounds is stabilized by the spin-density- 
wave instability due to Fermi-surface nesting or by 
anisotropic in-plane exchange interactions due to Sd or- 
bital orderingi^"— More importantly, spin fluctuations in 
the doped compounds, which are arguably a key to un- 
derstand the pairing mechanism, are largely unexplored. 
In particular, the question remains whether the normal- 
state spin fluctuations arc simply governed by the Fermi- 
surface topology, or other effects, such as orbital fluc- 
tuations. Indeed, it has been proposed that orbital-spin 
fluctuations in a multiband ground state could give rise to 



the superconducting pairing. ^^d^ Measurements of mag- 
netic fluctuations in the normal state should provide vital 
information to resolve these issues. 



Hence our goal is to investigate the change in the spin 
dynamics as a function of doping and elucidate the in- 
terconnection with the band structure. For the electron- 
doped compounds, recent ARPES and transport stud- 
ies clearly show the disappearance of the hole pockets 
around the antiferromagnetic zone center, i. e., a Lif- 
shitz transition.— This occurs for an electron doping xl 
of Ba(Fei„^CO:r)2As2 (0.15 < xl < 0.3 from ARPES, 
whereas ^ 0.1 from transport measurements) j^^— 
The doping dependence of the spin dynamics, however, 
has not been comprehensively studied on single-crystal 
samples by neutron scattering. Previous work focused 
on a powder sample^* or on the spin resonance in the op- 
timally doped and underdoped compounds. ^^^^"^ Here we 
investigate the spin dynamics in single-crystal, electron- 
doped Ba(Fei_:^Co^)2As2 for x = 0, 0.06, 0.14, and 0.24, 
ranging from the parent compound to the heavily over- 
doped regime, with a particular emphasis on magnetic 
fluctuations in the normal state. We find that the mag- 
netic fluctuations in the paramagnetic state of the parent 
compound are remarkably similar to those in the nor- 
mal state of the optimally doped compound. On the 
other hand, the spin dynamics in the overdoped regime 
is drastically different, and resembles that in the over- 
doped cuprates. As the cobalt content increases well into 
the heavily overdoped regime, magnetic scattering disap- 
pears, coinciding with the absence of superconductivity. 
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FIG. 1: (Color online) Constant-Q scans were measured on samples of four cobalt concentrations, (a) a; = at Q = (1,0, 1) 
(BT7), (b) X = 0.06 at Q = (1,0,1) (HB3), (c) x = 0.14 at g = (1,0,0) (GPTAS), and (d) x = 0.24 (HB3). Background 
is shown by open symbols. For x = 0, background was estimated from constant-energy scans at hcj — 3, 5, 8, and 15 meV. 
In (a), open red circles and open blue squares denote background calculated from constant-energy scans at 140 K and 250 K, 
respectively. For x = 0.06 and 0.14 background was measured away from the peak positions at Q = (1.2,0,1) and (3,0,0), 
respectively. For x = 0.24, both constant-Q and constant-energy scans do not show scattering intensity above background, 
which was measured at Q = (1.5,0, 1). The dotted lines in (a), (b), and (c) denote the fitted background 



II. EXPERIMENTAL DETAILS 

All single-crystal samples were grown from a self-flux 
using the Bridgman method described in Ref. [25l. The 
cobalt content was determined by energy dispersive x- 
ray analysis (EDX) using a scanning electron microscopy. 
Magnetic susceptibility of the x = 0.06 and x = 0.14 com- 
pounds exhibits superconducting transitions (onset) at 
26 K and 7 K, respectively. The transition temperatures 
place the former close to the optimally doped regime and 
the latter in the overdoped regime. For the x = 0.24 com- 
pound, the superconducting state is not observed down 
to 1.8 K. For each composition, up to four single crystals 
were co-aligned yielding a total mass of about 1 g. Inelas- 
tic neutron scattering measurements were performed on 
the triple-axis spectrometers BT7 at the NIST Center for 
Neutron Research, HB3 at the High Flux Isotope Reac- 
tor, Oak Ridge National Laboratory, and ISSP-GPTAS 
at the Japan Atomic Energy Agency. For clarity, all 
momentum transfers are labeled using the orthorhombic 
space group Fmmm, a low-temperature phase of the par- 



ent compound, even though the proper crystal structure 
is the tetragonal space group lA/mmm. The a; = 0, 0.06, 
and 0.24 samples were aligned in the hOl zone, while the 
X = 0.14 sample was aligned in the hkO zone. The final 
neutron energy was fixed at 14.7 meV. Pyrolytic graphite 
(PG) crystals were used to monochromate and analyze 
the incident and scattered beams using the 002 reflec- 
tion, respectively. Horizontal collimations of open— 50' — 
sample - 50'-open, 48' - 60' - sample - 80' - 120', and 
40' - 80' - sample - 80' - 80', were employed at BT7, 
HB3, and GP-TAS, respectively. PG fllters were placed 
in the scattered beam. The samples were cooled using a 
closed cycle ^He cryostat. 



III. MAGNETIC FLUCTUATIONS IN THE 
NORMAL STATE 

The scattering intensity can be written as SxiQ,!-^) = 
(n(a;,T) -I- l)xT'(q, w), where n(Ld,T) is the Bose factor 
and x" is the imaginary part of the spin susceptibil- 
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FIG. 2: Representative in-plane ((a)-(d)) and out-of-plane 
((e)-(h)) constant-energy scans at Tiu — 5 meV in the param- 
agnetic state of the parent compound. The solid lines denote 
the results of the global fit. 



ity. According to the theory of nearly antiferromagnetic 
metals2&, Xt(Q:'^) is given by 



Xo{T)T{T)huj 



{hujy -f r(T)2 . (1 + D^ql + F^{ql + ql)Y ' 

(1) 

where T is the damping constant, D and F represent 
the magnetic correlation lengths along the out-of-plane 
and in-plane directions, respectively, xo represents the 
isothermal susceptibility, and q=(qa,9b,'Zc) is a wave vec- 
tor away from an antiferromagnetic zone center. Back- 
ground was estimated from a series of constant-energy 
scans or a constant-Q scan taken away from the antifer- 
romagnetic wave vector and fit to a polynomial function 
of both momentum transfer and energy transfer. The co- 
efficients of this polynomial were initially assumed to be 
temperature dependent. 



A. Undoped t = 

For the parent compound {x — 0), our analysis of 
the background shows weak temperature dependence in 
the measuring temperature range. We, therefore, assume 
temperature-independent background in the fitting (see 
Figure 2). The background-subtracted scattering inten- 
sity St{Q,oj) was then converted to x" ■ Representa- 
tive constant-Q scans, constant-energy scans, St{Q,^), 
and x" measured in the paramagnetic state are shown 
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FIG. 3: Representative in-plane ((a) and (b)) and out-of- 
plane ((c) and (d)) constant-energy scans at hu — 9 meV in 
the normal state of the optimally doped compound. The solid 
lines denote the results of the global fit. 



in Figures (Ha), [51 and |lUa)-(b), respectively, x" 
fitted to the theory of nearly antiferromagnetic metals 
according to Eq. (1). In the fitting procedure, xo was 
constrained to obey the Curie law, since the data at high 
energy (huj ^15 meV), which are required to uniquely 
determine F and xo at high temperatures, are not avail- 
able. Error bars correspond to three times the standard 
deviation and in Figure [5] indicate large uncertainty at 
high temperatures, where F lies beyond the measuring 
energy range. We found that D and F do not change 
significantly within the measuring temperature range be- 
tween 140 K and 250 K, and the difference lies within 
the uncertainties. Therefore, we were unable to deter- 
mine their temperature-dependence and hence their val- 
ues were fixed at 2.6(5) A and 20(6) A, respectively. We 
note that the uncertainties were obtained from fitting 
the 140 K data and that they could become larger at 
high temperatures (see Figure [2|) . We have previously 
reported these anisotropic magnetic fiuctuations in the 
paramagnetic state of the parent compound at 145 K— . 
We have refitted those data using Eq. ([T]) (not shown), 
and have found that the lineshapes can be well described 
by Eq. ^ using the same parameters D and F, while F 
at 145 K was obtained from the linear relation mentioned 
below (also see Figure [5]). 

At each temperature, the solid lines in Figures [IJa), [5] 
and [4{a)-(b) denote the global fits to several constant- 
energy scans and the constant-Q scan, convoluted with 
the four dimensional resolution function. The resulting 
fit parameter F is linearly proportional to temperature, 
e.g., F(r) — a ■ T, where a = 0.16(6) meV/K as shown 
by the solid line in Figure O Note that F remains finite 
at the ordering temperature T/v = 136 K, and we do not 
observe any divergence of the correlation length at T/v. 
Therefore, unlike the parent compounds of the cuprates, 
the observed magnetically ordered state in BaFe2As2 is 
most likely not driven by the spin dynamics of the para- 
magnetic phase, but may be explained in light of the 
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FIG. 4: (Color online) The background subtracted scattering intensity and the imaginary part of the spin susceptibility x" 
at different cobalt concentrations, (a)-(b) a; = at Q = (1,0,1), (c)-(d) x = 0.06 at Q = (1,0,1), and (e)-(f) x = 0.14 at 
Q — (1,0,0). The dotted lines are guides to the eye. The solid lines denote the global fits to Eq. [1] convoluted with the 
resolution function. We note that x"{Qj<^) is identically zero for u — since it is required to be an odd function of to. 



interconnection between the lattice and magnetic inter- 
actions. 



B. Optimally doped x = 0.06 

For the optimally-doped compound {x = 0.06), the an- 
tiferromagnetic order is completely suppressed, and su- 
perconductivity emerges for T < Tc = 26 K. In the su- 
perconducting state, we observe the broad inelastic scat- 
tering centered at huj = 9.6 meV (Figures [TJb) andSfc)- 
(d)), in agreement with earlier reports j^i^ where the 
peak is attributed to the resonance mode. As tempera- 
ture increases above Tc, the inelastic peak is replaced by 
quasielastic magnetic fluctuations. Similar to the parent 
compound, the imaginary part of the spin susceptibility 
can be well described by the theory of nearly antiferro- 
magnetic metals [Eq. ([T])], which is shown by the solid 
lines in Figures [T]Jb), [3l and|4lc)-(d). The magnetic cor- 
relation length D (F) is equal to 2.4(6) A [19(3) A)] at 
30 K and 0.9(9) A [12(3) A] at 100 K. If measured at 
the same temperature, the correlation lengths measured 



in the optimally doped compound are shorter than those 
measured in the parent compound, which could be due to 
the change in the spin concentration upon doping. The 
in-plane magnetic correlation length is consistent with 
the earlier report on the x — 0.075 compound^*. We note 
that the out-of-plane magnetic correlation length was not 
measured in Ref. [28l . This anisotropic magnetic fluctu- 
ations observed in the paramagnetic state of the parent 
compound and in the normal state of the optimally doped 
compound are reminiscent of the anisotropic exchange 
interactions measured in the ordered state of the parent 
compound>22, 

We observe a marked similarity in the magnetic fluc- 
tuations in the normal state of the optimally doped com- 
pound and those in the paramagnetic state of the parent 
compound. More importantly, the temperature depen- 
dence of r follows the same linear relation as that ob- 
served in the parent compound (Figure [5]) . As a com- 
parison, Inosov et al.- - reported that the temperature 
dependence of F follows a similar linear form F(T) = 
a-{T+e) in the x = 0.075 compound (T^ = 25 K), where 
a = 0.14(4) meV/K and 8, the Curie- Weiss tempera- 
ture, is equal to 30(10) K. The fact that the imaginary 
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part of the spin susceptibility in the parent and optimally 
doped compounds can be well described by the theory of 
nearly antiferromagnetic metals suggest that the mag- 
netic fluctuations in the normal state of the parent and 
optimally doped compounds could have a common ori- 
gin, and are likely related to the presence of the quasi- 
two-dimensional hole and electron pockets observed by 
photoemission;^^ more evidences of this interconnection 
will be presented below. 



C. Overdoped x = 0.14 

In contrast to the parent and optimally doped com- 
pounds, the spin dynamics in the overdoped, supercon- 
ducting compounds [x — 0.14) shows the depletion of 
the scattering intensity in the low-energy region and gap- 
like excitations around 10 meV (Figures (He) and IH^e)- 
(f)). Its peak profile cannot be fit to either Gaussian 
or Lorentzian line shapes. Furthermore, the scatter- 
ing intensity exhibits weak temperature dependence, and 
the in-plane magnetic correlations are much shorter than 
those of the parent and optimally doped compounds as 
shown in the inset of Figure (He) . The imaginary part 
of the spin susceptibility [Figure Sff)] displays linear 
energy-dependence at low energy and a sharp drop at 
high energy. With the exception of the gap-like behav- 
ior, these magnetic excitations share many characteristics 
with those observed in the overdoped cupratesi^ 

The origin of the magnetic excitations in the over- 
doped, superconducting compound (a; = 0.14) is unclear 
at the present moment. We note that a; = 0.14 is a 
higher doping than the neck-collapsing Lifshitz transition 
(x ^ 0.1), and is indeed very close to the hole-pocket- 
vanishing Lifshitz point Thus, it is likely that the dom- 
inant part of the hole band lies below the Fermi level. Re- 
cent theoretical calculations show that in such a case the 
imaginary part of the spin susceptibility is strongly sup- 
pressed giving rise to the pseudogap behavior'^", and such 
a pseudogap was observed in NMR measurements in the 
electron overdoped regimes of LaFeAsOi_a;Fa; (Ref. Isil ) 
and Ba(Fei_a;Coa;)2As2i^ Even though the gap-like be- 
havior observed in our neutron scattering suggests that 
the majority of the hole Fermi surface already diminishes 
at X = 0.14, superconductivity with lower Tc = 7 K is 
still observed. This superconductivity may then be of the 
nodal type, where pairing is formed between electrons on 
the same electron Fermi surface around the zone corners, 
as has been theoretically proposed as one alternative to 
the s± mechanismiSS, A recent heat transport experiment 
indeed indicates that the superconducting gap shows a 
tendency to be strongly anisotropic as the electron dop- 
ing increases in Ba(Fei_2:Co2:)2As2r^ although the mea- 
surements were made only up to a; = 0.114 so a direct 
comparison with our result (x = 0.14) is not possible at 
present. 




T[K] 

FIG. 5: (Color online) The damping constant F as a function 
of temperature. Red circles and blue squares represent the 
damping constants of the parent and optimally doped com- 
pounds, respectively. The dotted line shows the linear relation 
measured on the optimally doped compound {x = 0.075)2^. 
The solid line denotes our result (see the text). 



D. Heavily overdoped x = 0.24 

For the heavily overdoped, non-superconducting com- 
pound [x = 0.24), Figure [TJd) shows the suppression 
of the magnetic scattering. NMR measurements on the 
X = 0.26 compound"^^ and neutron scattering measure- 
ments on electron-doped LaFeAsOi-jjOa; in the heavily 
overdoped regimeii^ reveal the suppression of the spin 
fluctuations consistent with our result. In this regime, 
both photoemission measurements and first-principles 
calculations point to the disappearance of a hole Fermi- 
surface pocket.— Our result, therefore, further sug- 
gests the correlation between the electronic band struc- 
ture and magnetism, and supports the scenario that the 
magnetic fluctuations in the underdoped and optimally 
doped regimes, which serve as a precursor to supercon- 
ductivity, originate from quasiparticle scattering across 
the electron and hole pocketsi^ 



IV. CONCLUSION 

We have studied the spin dynamics in electron-doped 
Ba(Fei„a;Co2;)2As2 at four cobalt concentrations. We ob- 
serve a striking similarity between the magnetic fluctua- 
tions in the paramagnetic state of the parent compound 
and those in the normal state of the optimally doped 
compound, and the suppression of the magnetic signal in 
the heavily overdoped regime, in which superconductiv- 
ity disappears. These two results suggest that magnetism 
and superconductivity are strongly correlated. On the 
other hand, magnetic excitations in the vicinity of the 
hole-pocket-vanishing Lifshitz point are markedly differ- 
ent from those in the underdoped regime, with the emer- 
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gence of a spin gap and much weaker temperature de- 
pendence of x"- These changes in the spin dynamics at 
different doping levels likely reflect changes in the elec- 
tronic band structure. Further experimental and theo- 
retical work is desirable to examine the interconnection 
between the band structure, spin fluctuations, and su- 
perconducting pairing mechanism. 
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